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Abstract 



Perfect fluid spacetimes admitting a kinematic self-similarity of infinite type are 
investigated. In the case of plane, spherically or hyperbolically symmetric space-times 
the field equations reduce to a system of autonomous ordinary differential equations. 
The qualitative properties of solutions of this system of equations, and in particular 
their asymptotic behavior, are studied. Special cases, including some of the invariant 
sets and the geodesic case, are examined in detail and the exact solutions are provided. 
The class of solutions exhibiting physical self-similarity are found to play an important 
role in describing the asymptotic behavior of the infinite kinematic self-similar models. 



PACS numbers: 04.20.Ha, 04.20. Jb, 04.40.Nr, 98.80.Hw 



1 Introduction 



Since the pioneering work of Sedov [|IJ, the study of self-similar systems has played an im- 
portant role in an extensive range of physical phenomena in the classical (Newtonian) theory 
of continuous media, giving rise to many very interesting results with useful experimental 
and astrophysical applications. A characteristic of self-similar solutions is that, by a suit- 
able transformation of coordinates, the number of independent variables can be reduced by 
one, thus allowing a reduction of the field equations (eg., in some cases partial differential 
equations can be reduced to ordinary DEs). Such solutions are of physical relevance since 
they are often singled out from a complicated set of initial conditions; for instance, in an 
explosion in a homogeneous background solutions asymptote to self-similar solutions. 

In the context of general relativity, the concept of self-similarity is also largely doc- 
umented in the literature, beginning with the pioneering paper by Cahill and Taub ||, and 
followed by important work by Eardley [|], [5|. Spherically symmetric homothetic solutions 
were studied, which proved to be especially useful in the cosmological context. More re- 
cently Carter and Henriksen || |7| have introduced the concept of kinematic self-similarity, 
which is a generalization of the homothetic case. 

The existence of self-similar solutions of the first kind (homothetic solutions) is re- 
lated to the conservation laws and to the invariance of the problem with respect to the group 
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of similarity transformations of quantities with independent dimensions. In this case a cer- 
tain regularity of the limiting process in passing from the original non-self-similar regime 
to the self-similar regime is implicitly assumed. However, in general such a passage to this 
limit need not be regular, whence the expressions for the self-similar variables are not deter- 
mined from dimensional analysis of the problem alone. Solutions are then called self-similar 
solutions of the second kind. Kinematic self-similarity is an example of this more general 
similarity. Characteristic of these solutions is that they contain dimensional constants that 
are not determined from the conservation laws (but can be found by matching the self-similar 
solutions with the non-self-similar solutions whose asymptotes they represent) 0. 

In the study of relativistic dynamics, there is an important distinction which must 
be made. The existence of a symmetry for the geometry (i.e. the metric) does not neces- 
sarily imply the existence of a symmetry for the matter functions (in particular, the energy 
density and pressure, when considering a perfect fluid). For that reason it is important to 
distinguish between the ideas of "physical" self-similarity and "geometrical" self-similarity. 
The definitions have been given elsewhere ||, and it has been shown that in the case of 
finite kinematic self-similarity the subclass of solutions exhibiting physical self-similarity 
have an important role to play in the examination of the full dynamics |H| |H| . Similar 
investigations will also be important in the study of infinite kinematic self-similarity, and 
this is what we shall consider here. 



In Benoit and Coley |ID] it was shown that in the case of finite kinematic self- 
similarity all solutions which exhibit physical self-similarity asymptote (in past and future) 
to solutions which exhibit a homothety. In fact, this result was extended in Benoit, |T]J to 
consider all solutions with finite kinematic self-similarity (i.e., not simply those exhibiting 
physical self-similarity). 

The definitions given previously || [7], [TU], |TTJ for the relativistic kinematic self- 
similarity show the dependence on a parameter, commonly denoted by "a". In the work 
of Benoit and Coley [pi], |TTJ the value of this parameter is assumed to be an arbitrary fi- 
nite value. One special case that was not considered in the previous papers is the case in 
which a takes on an 'infinite value'. This case corresponds to the generalization of rigid 
transformations in general relativity. What we attempt in this paper is to study those space- 
times containing a non-null 2-space of constant curvature when they admit a kinematic self 
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similarity of infinite type, thus complementing the previous work. We shall place special 
emphasis on models which can be interpreted as perfect fluid solutions of Einstein's field 
equations (EFE). In this case the governing system of differential equations reduces to a 
system of autonomous ordinary differential equations and we shall analyze the qualitative 
behavior of these models. Exact solutions are obtained in some special cases, particularly 
those which are of importance in the asymptotic analysis. 

The paper is organized as follows: Section 2 contains a brief description of kinematic 
self-similarity, deriving the form of the kinematic self-similar vector field and the self-similar 
equations for space-times admitting a three-dimensional multiply transitive group of isome- 
tries. Section 3 contains the details of the reduction of the EFEs when there exists a proper 
kinematic self-similar vector that commutes with all of the Killing vectors. The equations are 
considered in the different cases characterized by the orientation of the fluid flow. Section 4 
examines the nature of solutions to these equations through the use of qualitative methods. 
Section 5 provides the physical asymptotic solutions. Special cases are then studied in more 
detail in sections 6 and 7. 



2 Kinematic self-similarity and perfect fluids 

A vector field £ is called a kinematic self-similar vector (KSS) if it satisfies the conditions 
i 

C£U a = au a and C^h ab = 26h ab , (2.1) 

where a and 5 are constants, £ stands for the Lie derivative operator, u a is the four-velocity 
of the fluid and h ab = g a b + u a u b is the projection tensor which represents the projection 
of the metric into the 3-spaces orthogonal to u a . Evidently, in the case a = 5 it follows 
that £ is a homothetic vector (HV) corresponding to a self-similarity of the first kind, and 
ifa = 5 = 0,£isa Killing vector (KV). 

The similarity transformations are characterized by the scale-independent ratio, a/S, 
which is referred to as the similarity index. This index is finite except in the case of rigid 
transformations characterized by 5 = 0. In this case the self-similarity is referred to as 
'infinite' type. Further information regarding KSS and their properties can be found in 
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This paper focuses on kinematic self-similar models exhibiting a three-dimensional 
multiply transitive group of isometries, G3. Since we consider only perfect fluid models, the 
G3 necessarily acts on space-like orbits S2. The solutions then correspond to spherical, plane 
and hyperbolic symmetric space-times, and the line element of the metric can be written in 
comoving coordinates as 

dg 2 = _ e 2*(^)^2 + e 2*(t,r) dr 2 + 5 2 ( ^ r)( ^2 + E( ^ A . ) 2 #2) ? ^ 

where 

{sin# k = +l 
9 k = (2.3) 
sinh 9 k = — 1 . 

The four-velocity vector is then given by 

« a = (- c »(*.0 0,0,0) • (2.4) 

The Killing vectors (KVs) for the space described by the metric ( j2.2p are 

£' 

771 = sinrftde + —cos 

Zj 

£' 

i]2 = cos(f>dg — —sin 

Zj 

where a dash denotes a derivative with respect to 9. These KVs satisfy the following 
commutation relations: 

fax, 7] 2 ] = kri 3 , [r}2, 7] 3 ] = 771 , [773, 771] = T] 2 . 

If we assume the existence of a proper KSS, £, then from the Jacobi identities and 
the fact that the Lie bracket of a proper KSS and a KV is a KV, the following algebraic 
structures are possible 

(/) [&tfc] = 0, i = 1,2,3, fc = 0,±l, 
€ = ^(t,r)^ + f(t,r)9 r , 

(//) = 771 , [tm] = V2, [£,m] = o, k = o, 

£ = e(t,r)d t + e(t,r)d r -6d e . 
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Focusing attention now on equations fl2.1|) for the KSS and the particular metric (|2.2|) , it is 
easy to show that the KSS takes the form 

£ = ?(t)dt + e{r)d r + m0d , (2.5) 

where m is a real number. 

In the case in which £ commutes with all of the KVs, i.e., m = 0, the self-similar 
equations (|2.1|) reduce to 

ut,^ + u t,rC + u tC\t = au t (2.6) 

I fa/ yip I ^ifa/yp^ — 2(5/iy^* ^2.7^ 

V*6* + Vrf = 25/i ee , (2.8) 

where a comma indicates partial derivative. In this case, the three metric forms (plane, 
spherical and hyperbolic) can be studied together. The form of the metric functions are 
similar and the EFE's reduce to single system of ODE's. The second algebraic structure, 
(II), is only possible in the plane symmetric case, in which case the form of metric functions 
and the governing equations are different from those studied here. This last case will be 
studied elsewhere. 

In the infinite case, 5 = and £ can be normalized so that the constant a can be set 
to unity, which we shall do hereafter. 



3 Reduction of EFEs 

When attention is restricted to the case in which the KSS, commutes with all the KVs, 
(i.e., there exits a proper KSS, orthogonal to all the Killing vectors) and in which the 
KSS is of infinite type, three different cases arise. The three cases are dependent on the 
orientation of the fluid flow u relative to the KSS; i.e., fluid flow parallel to fluid flow 
orthogonal to and the most general 'tilted' case. Each case will now be discussed, with 
the focus and the detailed analysis made in the general 'tilted' case. 
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3.1 Fluid flow parallel to £ 



In this case £ takes the form £ = £*(i)9t and without loss of generality we can choose it to 
be £ = tdf. The metric can then be written as 

d s 2 = _ e 2Hr) dt 2 + dr 2 + S 2 (r)(d6 2 + E(0, &) 2 d0 2 ) . (3.1) 
The EFEs for the perfect fluid become 

^S 2 = k - S 2 r - 2SS t „ (3.2) 
pS 2 = —k + S 2 r + 2SS >r <b jr (3.3) 

= ^ + + + (3-4) 

where k — 1,0, —1 for spherical, plane and hyperbolic symmetry, respectively. (Note that d t 
is a KV). Equation ( |3.1| ) represents a static space-time. The function S )T vanishing implies 
fM+p — 0, and therefore realistic perfect fluid solutions are excluded. Apart from this case, 
any functions S and $ satisfying (|3.4j) represent kinematic self-similar solutions. 



3.2 Fluid flow orthogonal to £ 

In this case £ takes on the form £ = £ r (r)d r , and without loss of generality we can choose 
it to be £ = d r . The metric can then be written as 

ds 2 = -e 2r dt 2 + e 2 * ( *W 2 + S 2 (t)(d6 2 + E(0, kfd<p 2 ) . (3.5) 

The field equations for a perfect fluid are 

= S >t (3.6) 
= 1 + ke™S- 2 - e 2 *~ 2r (^ 2 t + * >tt ) (3.7) 
^ = fcS" 2 , p= -A;^ 2 . (3.8) 

This case is again empty of perfect fluid solutions with + 
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3.3 General 'tilted' case 



The general 'tilted' case occurs when the four-velocity is neither parallel nor orthogonal to 
the self-similar vector field. In this case one can choose coordinates so that the KSS takes 
the form: 

£ = td t + rd r . (3.9) 

In such coordinates, and solving equations ( |2.6|) -( [2~8"|) , it is easy to show that the metric can 
be given by 

ds 2 = -e 2 *dt 2 + —dr 2 + S 2 {d6 2 + E(0, kfdcf 2 ) , (3.10) 

where $, \I> and S are functions depending only on the self-similar coordinate 

t 



(3.11) 



The field equations for a perfect fluid source are 



= SH + ffr-Q'S't-V'S't 
= ^(O + E^), 



(3.12) 
(3.13) 



where a dash denotes derivative with respect to £ and 



Si 



JS') 2 + -S 2 & - S 2 W + S 2 ^" + k^-e 2 * + S 2 (<!>') 2 
t ' 



r 
t 2 





2 = e 



2*-2$ 



SS' - S 2 {^') 2 + S 2 &V' - S V + {S 



-f/\2 



(3.14) 
(3.15) 



The only possible solutions to equation ( |3.13j ) must necessarily satisfy Ei = E 2 = 0. 



Assuming \x + p ^ 0, we have that 5" cannot vanish. Then defining z = ln(£) and 
/ = df/dz, equations (|3.12|) and ( |3.13| ) can be rewritten as 









S 

s 



. $ _ q, 

I)' 



„2* 



+ k— + $ 2 



~s + [s. 



\Ir + - * + \& . 



(3.16) 
(3.17) 

(3.18) 
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Applying now the definitions (to be consistent with the notation in [ffO ) 

y = S/S , u = 4f , v = ^> , w = -ke 2 *S~ 2 , (3.19) 
equations (|3.16|) - (|3.18|) reduce to a 4-dimensional autonomous system of ODEs 



y = y( u + V -y) (3.20) 

u = —y + y 2 — u 2 + u + uv (3-21) 
v = w + y 2 + uv - v 2 (3.22) 
w = 2w{u - y) . (3.23) 

The matter quantities are given by 

H = e- 2 *[2yu + y 2 }r 2 - e' 2 ^ [w + 2yv + y 2 ] (3.24) 
p = e -^[-2yu + 2y-y 2 ]r 2 + e~ 2 * [w + 2yv + y 2 }. (3.25) 

Note that the density and pressure can be split as fi = fii + fi 2 and p = p\ + p 2 where fii = 
£ii((,)t~ 2 , pi = Pi(Qt~ 2 and — p 2 — ^2 — /^(O- Each component of the density and pressure 
then exhibits self-similarity in that C^ji\ = — 2/xi, Cgpi = —2p\ and C^fi 2 = C^p 2 = 0. 

We note the following special cases which are evident from equations ( |3.24| ) / ( ft.25| ): 

1. In the particular case w + 2yv + y 2 = (i.e., ji 2 = p 2 = 0) the fluid is said to be 
'physically' self-similar ||. 

2. The case y = is also 'physically' self-similar; and since /i + p = solutions in this 
case correspond to a cosmological constant solution. 

3. The case y = w = gives rise to vacuum solutions. 

4. Perfect fluid solutions (with /i + p 7^ 0) will exhibit a barotropic equation of state 
(p = p(n)) if and only if 

w + 2yv + y 2 = Coe 2 * and 2u + y = C\ , (3.26) 

where c and C\ are constants. 
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5. If we are to demand that the solutions satisfy the weak and dominant energy conditions 
(i.e. — p < p < p) over the entire manifold the following inequalities serve as necessary 
conditions 

y>o , 

y(2u + y-l)>0 , 

w + 2yv + y 2 < . (3.27) 

We therefore note that by demanding the energy conditions be satisfied throughout 
the evolutions of these models, the possible asymptotic behaviors are greatly reduced. 

Each of these cases will be important in the analysis of the equations, which follows in the 
next sections. 



4 Qualitative analysis 



The system given by equations (|3.20|) - (|3.23|) is an autonomous system of first order ODEs. 
As such, the asymptotic behavior of the system can be determined by studying the quali- 
tative dynamics. 



The full system of equations, ( |3.20|) - (|3.23| ), describing all possible solutions, exhibits 
a number of invariant sets, including the planes 



Ji : w = , J 2 : y = u , J 3 : y = 



as well as the surfaces 



Z4 : w + 2yv + y 2 = , I5 : w — yv + y 2 = . 
To allow for the simplification of the analysis we make the following change of variables: 

xi = y, x 2 = u- y, x 3 = v, x 4 = w . 



In these coordinates the equations (|3.20 ) - (|3.23| ) become 



Xl 



Xi{X 2 + x 3 ) 



(4.1) 
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x 2 = x 2 (l + x 3 - x 2 - 3xi) , (4.2) 
x 3 = x 4 + x\ + X1X1 + x 2 x 3 - xl , (4.3) 
x 4 = 2x 4 x 2 . (4.4) 

The finite singular points can then be located (note, this system is not bounded). 
They are summarized in Table |I|. There are three distinct hyperbolic singular points and 
two sets of non-isolated singular points, each of which have zero eigenvalues in the direction 
tangent to the curve and non-zero eigenvalues in all other directions (i.e., they are normally 
hyperbolic). The finite singular points can be classified by the eigenvalues of the Jacobian 
for the vector field. This classification is given in Table [I], and will be discussed in the 
sections to follow. 

We can also consider the singular points located at infinity. To do this we employ a 
Poincare transformation using the variables: 

Xi = Xi9, X 2 = x 2 9, X 3 = X36 1 , X 4 = x^6 , 
9 = (1 + xj + x\ + x\ + xfj- 1/2 . 



In this case the equations (|4.1| )- (|4.4| ) become 



X[ = X 1 (X 1 + X 3 -K)-X 1 (Xi + X 3 X i )9, (4.5) 

X' 2 = X 2 (X 3 - X 2 - 3X 1 - K) -X 2 (Xl + X 3 X 4 -X 2 )9, (4.6) 

X' 3 = X 2 1 +X 3 (X 1 +X 2 -X 3 -K)-X 3 (X 2 2 + X 3 X 4 -X 4 )9 } (4.7) 

X> 4 = X 4 (2X 2 -K)-X 4 (X 2 + X 3 X 4 )9, (4.8) 



where 



K = X\X 2 + 1X\X 3 — ?)X\X\ + X\X^ — X\ 

+ X 2 2 X 3 + X 2 Xl - Xl + 2X 2 X 4 . (4.9) 

The singular points located on the invariant boundary 9 = [the location of the 
infinite singular points for equations ( |4.1|) -( f4~4D l can be classified by examining the dynamics 
restricted to this invariant surface. The location of the singular points and their classification 



are identical to that of the finite case |TTJ, and are given in Table 
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Returning now to the system ( fO) -(O), we see that the invariant hyperplanes X4 = 
and x 2 = divide the phase space into four additional invariant sets: 



51 = {(x 1 ,x 2 ,x 3 ,x 4 )\x 2 > 0,x 4 > 0} , (4.10) 

5 2 = x 2 ,x 3 , x A ) \x 2 > 0,x 4 < 0}, (4.11) 

53 = {(x 1 ,x 2 ,x 3 ,x 4 )\x 2 < 0,x 4 > 0} , (4.12) 
5" 4 = {(xi, x 2 , x 3 , Xi)\x 2 < 0, X4 < 0} . (4-13) 

In each of these invariant sets, the function X4 (curvature) is monotonic. As a result all 
stable asymptotic behavior is necessarily located on one of the invariant sets x 2 = or 
£4 = (or at £4 = ±00). Each of these cases will be studied separately. In all cases we note 
that the classification of the singular points (both finite and infinite) can be determined by 
considering the points listed in Tables [l] and § restricted to the invariant set being considered. 



4.1 Subcase: X2 = 

We first consider the hyperplane x 2 = 0. In this case the system of equations ([4. 1|) - ( |4.4j ) 
becomes: 

x\ = x x x 3 , (4.14) 
^3 = ^0 + ^1 + ^1^3-^3, (4-15) 
x\ = 0; £4 = w = const . (4.16) 

This system is a two-dimensional dynamical system in the variables x\ and x 3 with parameter 
Wq. The finite singular points are located (where they exist) at: 

L 1± = (0,±V^), (4.17) 
L 2± = (±7=^,0). (4.18) 

Each of these points is the intersection of the fixed curves (L± and L 2 respectively) with the 
plane under consideration; i.e., x 2 = and x 4 = w . 

We note here that the value wq = is a bifurcation. We shall first consider the 
dynamics of the solutions when wq < and wq > 0, considering the dynamics at the 



11 



bifurcation point after. We can see from Table [T] that when wq < or u> > 0, the points 
Li± and L 2 ± have both positive and negative eigenvalues when restricted to this case. 
Therefore, each is a two-dimensional saddle point. Further, from Table [2] we see that the 
infinite singular points are B±, C±, D±, and E±. Restricted to this invariant set we find 
that B + , D_ and E_ are sources; whereas B_, D + and E + are sinks. 



We now turn our attention to the dynamics of equations ( }4.14|) -( f4.15|) at the bifur- 



cation value of Wq = 0. In this case we see that there is only one finite singular point, which 
is located at the origin, (21,23) = (0,0), i.e., at the intersection of the two fixed curves L\ 
and L 2 . This singular point is non- hyperbolic in nature, and as such its local properties can 
not be determined by examining the eigenvalues of the corresponding Jacobian matrix. In 
this case, however, there are three invariant lines: x\ — 0, x\ — —2x 3 and x\ = x 3 . The 
dynamics on each of these lines can be determined as follows: 

(i) on Ji : x\ = 0: x 3 = —x\ < 0. 
(ii) on J 2 : X\ = —2x 3 : x 3 = x\ > 0. 
(iii) on J 3 : x\ = x 3 : x 3 = x\ > 0. 

Each of these three invariant lines then divide the 2-dimensional phase space into 6 additional 
invariant regions: 



Ji 


= {(xi, 


X3) 


\xi > 0, x 3 > Xi}: 


x\ 


> 


h 


= {(xi, 


X3) 


\x\ < 0, x 3 > —X\/2}: 


x\ 


< 


■h 


= {(xi, 


X3) 


|xi > 0,x 3 < -xi/2}: 


x\ 


> 


■h 


= {(xi, 


X3) 


\xi < 0, x 3 < Xi}: 


x\ 


< 


h 


= {(xi, 


x 3 ) 


\x\ > 0, —Xi/2 < x 3 < X\}\ 


x 3 


> 


■h 


= {(xi, 


x 3 ) 


\x\ < 0,Xi < x 3 < — xi/2]: 


x 3 


> 



The result is that the point (0, 0) is a saddle. The asymptotic analysis is then completed 
by considering the singular points on the infinite boundary. As the quadratic portion of the 
vector field is unchanged by the differing values of the bifurcation parameter, the infinite 
singular points and the corresponding analysis is identical to that when wq ^ 0. 

A bifurcation diagram, including all the phase portraits for each range of the param- 
eter w is given in Figure |1|. As can be seen by these phase portraits all generic asymptotic 
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behavior (to the past and the future) is located on the infinite boundary. The exact solu- 
tions for each of these singular points (which are asymptotic states to past or future or are 
intermediate states) will be examined in section 5. 



4.2 Subcase: £4 = - plane symmetry 

The invariant set £4 = contains a subset of the asymptotic solutions for the system ( |4.1| )- 
( |4.4|) . As can be seen from equations ( |2.2|) and ( |3.19| ), solutions which have w identically 
zero comprise the set of plane symmetric solutions. 

In this case, the system of ODEs (|4.1|)-( fO| ) reduces to: 

x\ = xi(x 2 + £3) , (4.19) 
x 2 = x 2 (l + £3 - 3xi - x 2 ) , (4.20) 
£3 = x\ + x 3 (x 1 + x 2 - x 3 ) . (4.21) 



The co-ordinate planes x\ = and x 2 = are each invariant sets for this system, as are the 
sets Xi + 2x 3 = and X\ = x 3 . 

The finite singular points in this case are given by Qi, Q 2 , Q 3 and L\ = L 2 . The local 
dynamics of each is determined by considering the sign of the eigenvalues of the Jacobian (see 
Table 0) restricted to this set X4 = 0; i.e, those eigenvalues whose associated eigenvectors 
have the form (ci, c 2 , c 3 , 0). The points Q1-Q3 are saddles in this three-dimensional set. The 
point L = Li = L 2 is non-hyperbolic. Center manifold theory |13| allows the point to be 
analyzed. The many invariant sets which include this point greatly simplify the analysis, 
and it is a straightforward matter to show that in the two dimensions which define the 
coordinate plane x 2 = the point is a saddle and in the third direction it is a saddle-node. 
The infinite singular points (not including C±) are given in Table |2]. The dynamics on the 
infinite boundary is represented by Figures |2] and ^ (see [jTT| for details). 

Before considering the global dynamics in this three-dimensional system, we shall 
consider the dynamics as restricted to the invariant planes. Each of these planes will divide 
the phase space further, allowing for a simplification in the analysis when considering the 
entire space. Note that the x 2 = invariant set has been completely analyzed in the previous 
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section. The dynamics are represented by the case Wq = in Figure [T]. Therefore, we need 
only consider the planes x% — 0, X\ + 2x 3 = and X\ = x 3 . 

Invariant Set: x\ = 

In the invariant set X\ — 0, the system ( |4.19| )-( f4.21| ) reduces to: 



x 2 = x 2 (l + x 3 -x 2 ) , (4.22) 
x 3 = x 3 {x 2 - x 3 ) , (4.23) 

and represents the vacuum solutions in the full 4-dimensional system. This system gives rise 
to dynamics in the x 2 - x 3 plane. The finite singular points are given by L\ = L 2 = (0, 0) and 
Qi = (1, 0). Local analysis shows that the point (1, 0) is a saddle point and the point (0, 0) 
is non-hyperbolic, saddle-node in nature (determined through the use of center manifold 
theory). Therefore, no stable asymptotic behavior is located in the finite part of the phase 
space and all asymptotically stable solutions in this subcase are located on the infinite 
boundary. The complete phase portrait, as compactified by the Poincare transformation, is 
given in Figure |j. 

Invariant Set: x\ + 2x 3 = 

In the invariant set Xi + 2x 3 = 0, the system ( |4.19| )-( |4T21| ) reduces to: 



x 2 = £2(1 + 3x3-22), ( 4 - 24 ) 
x 3 = x 3 (x 3 + x 2 ). (4.25) 

As a result the dynamics is located in a two-dimensional plane. The finite singular points are 
given by Q\ = (1, 0), Q 3 = (1/8, —1/8) and L\ = L 2 = (0, 0). Local analysis determines that 
the point (1, 0) is a saddle, (1/8, —1/8) a spiraling sink and (0, 0) a saddle-node (determined 
through the use of center manifold theory) . The phase portrait for this compactified 
by the Poincare transformation, is given in Figure |5|. In this case the fluid is also physically 
self-similar. 

Invariant Set: X\ = x 3 

In the invariant set X\ = x 3 , the system ( |4.19| )-( |4~2""[| ) reduces to: 

x\ = xi(x! + x 2 ), (4.26) 
x 2 = x 2 {l-x 2 -2x x ). (4.27) 
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As a result the dynamics are located in a two-dimensional plane. The finite singular points 
are given by L\ = L 2 = (0, 0), Qi = (0, 1) and Q 2 = (1, —1)- Local analysis determines that 
the points (1, 0) and (1, —1) are saddles whereas (0, 0) is a saddle-node (determined through 
the use of center manifold theory). The phase portrait for this CclS6 ; clS compactified by the 
Poincare transformation, is given in Figure || 

Global dynamics 

The global dynamics can be determined through an investigation of the direction fields 



making use of the monotonicity principle |L4]]. To simplify the global analysis, consider 
the full three-dimensional phase space divided into 16 invariant regions and labelled Ui so 
that £/i3 corresponds to the set x 2 = 0, Uu the set x\ = 0, Ui$ the set x\ + 2x 3 = and 
£/i6 the set x\ = 23. In each of the remaining 12 regions of space a monotonic function 
has been identified. These regions, and their corresponding monotonic functions, are given 
in Table |3|. Note that the totality of the sets Ui, i = 1..16 provides a decomposition of 
the complete phase space. As such, since each region is invariant under the system ( |4.19| )- 
( |4.21| ), the existence of strictly monotonic functions in the regions U\- U14 ensures that the 
only possible asymptotic solutions are located on the boundaries (either finite or infinite). 
The finite boundaries are the sets Ui 3 — U%q, or subsets thereof. As a result the global 
dynamics has been completely determined by the previous investigations. The only possible 
asymptotic states are, therefore, the singular points located at finite and infinite values. 
Furthermore, in the full four dimensional space the only possible asymptotic states are 
those singular points which are sinks or sources; namely the sinks A + , B_ and D + and the 
sources A_, B + , and _D_. 



4.3 The case £4 = ±00 

^From table ^] we see that the only solutions characterized by 24 — > ±00 correspond to the 
points C±, found by compactifying the phase space using a Poincare transformation. This 
point is 'non-hyperbolic' in all four directions. To determine the exact nature of the local 
behavior of these points we consider the system ( f4.1|) -( [4~4|) under the following change of 
coordinates: 

y 1 = Ei f Y 2 = ^, F 3 = -, and Y 4 = -, (4.28) 
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and a "time" variable that is defined by /' = Y^f . The singular points of interest (namely 
C±) are now located at the origin of the new coordinate system. In these new coordinates, 
the system fl4.1|)- ([Op becomes: 

Y{ = Y 1 (Y 3 -Y 2 ) (4.29) 

y 2 ' = Y 2 (Y^ + Y 3 -3Y 1 -3Y 2 ) (4.30) 

Yi = Y, + Y 1 2 + Y 1 Y 3 -Y 2 Y 3 ~Y 3 2 (4.31) 

Y; = -2Y 2 Y 4 . (4.32) 

The singular points of interest (namely C±) are now located at (0,0,0,0). There are two 
invariant lines, namely Y\ = Y 3 = Y^ = and Y\ — Y 2 — Y± — 0. Each of these lines 
corresponds to an eigenvector of the flow for the system ( |4.29p -( fl.32| ), and on each of these 
lines the flow is monotonic decreasing. The dynamics in a third direction can then be 
determined by considering the two-dimensional set Y 2 = Y4 = 0, i.e.: 

Y{ = F1F3 (4.33) 
Yi = Y 2 -Y 3 2 + Y 1 Y 3 . (4.34) 

In this case Y\ = Y 3 and Y\ = —2Y 3 are invariant sets (and, again, eigenvectors of the flow). 
On each of these sets the derivatives are strictly positive or strictly negative, indicating that 
this point is a saddle-node. 

Therefore in three of the four directions (of the full phase space) through the singular 
points C± the derivative does not change sign, and in the fourth direction there is no 
motion (as this direction is normal to the sheets of invariant planes described in the previous 
section). This point is, therefore, a higher- dimensional saddle-node. 



5 Description of asymptotic solutions 

In the qualitative analysis of the previous section, the asymptotic states of the governing sys- 
tem were described as singular points of the autonomous system of ODEs. The existence of 
other types of stable structures was ruled out by the existence of monotonic functions. While 
some of the singular points being described are not structurally stable in all 4-dimensions, 
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there are invariant regions in which they do act as attractors (either to the past or the 
future); in addition, these points act as intermediate attractors (repellors) for large classes 
of solutions. Each of the physical solutions described by these singular points will now be 
given, restricting attention to those solutions which satisfy the weak and dominant energy 
conditions. In the case of infinite kinematic self-similarity, two of the boundaries of the 
regions which satisfy the energy conditions are, in fact, invariant sets; therefore we need 
only consider the solutions which lie in the regions x\ > and £4 + 2x1X3 + x\ < 0. For the 
invariant set X\ = 0, one has fi+p = corresponding to a cosmological constant or vacuum 
if, in addition, x 4 = 0. 



5.1 Finite singular point asymptotic states 

Qi = (0,1,0,0) is a vacuum solution corresponding to Minkowski space-time. Q 2 = 
(1,-1,1,0) does not satisfy the energy conditions since \x — p < 0. L\ — (0, 0, /?, /3 2 ), if 
(3 7^ 0, then \i + p = and /i < violating again the energy conditions. The remaining 
singular points are: 

. Q 3 = (1/4, 1/8, -1/8,0) 

In this case the metric is plane symmetric 

/ r \l/4 h2 /f\3/4 //\V2 

d* 2 = -Q) de + b -(^-) rfr 2 + (l) (d9* + 9W), (5.1) 
where b is a constant. The energy density and pressure are given by 

» = P=^. -) ■ (5-2) 



At 2 Vr, 

The fluid is physically self-similar and it represents the only stiff-matter solution (i.e., 
fi = p) in the plane symmetric case. 

L 2 = (a,0,0,-a 2 ) 

The case a = (i.e., the intersection of L x and L 2 ) corresponds to Minkowski space- 
time. Otherwise, the metric is spherically symmetric 

ds 2 = -dt 2 + ^ (£) 2U dr 2 + s 2 ^ (d9 2 + sm W) , (5.3) 
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where s and a are constants. The energy density and pressure are 

3a 2 -3<7 2 + 2ct 
» = ~jT > P = ^ • ( 5 - 4 ) 

In this case the fluid is physically self-similar and satisfies the energy conditions for 
a > 1/3. The case o = 1/3 corresponds to stiff- matter and a = 2/3 to dust. 



5.2 Infinite singular point asymptotic states 

The infinite singular points are displayed in Table |]. A±, B± and F± correspond to vacuum 
solutions. They asymptote respectively to the following line elements: 

ds\A) = -dt 2 + -J^ (t ^ ' j2 ^ 2 ' ^ ' ° 2 ^ 2 



+ 62 










In 




+ c 











dr z + d6 z + 



ds 2 (B) 



ds 2 {F) 



exp 



In ( -\ + c 

r J 

tV 2 ' 
<• I - 

r 



. for ln(-)+r — 0. (5.5) 



for In ( - ) : r — . (5.6) 



,2 fe 2 * 

ctt H — - exp 



,i<Y /2 

r J 



dr 2 +d6 2 +6 



2ji2 



for >■ oo , (5.7) 

r 



where 6 and c are integration constants. 



The singular points C+, -D-t, G + , and if-t do not satisfy energy conditions. Phys- 
ical solutions (i.e., perfect fluid solutions satisfying the energy conditions) do not asymptote 
to these infinite singular points since they lie in different invariant regions of the phase 
space. The remaining singular points are: 



C-\ The metric is given by 



ds 2 = -dt 2 + ^-dr 2 + S 2 (d6 2 + stn 2 ° < 2 



where 



In ( - ) + c 

r 



2/3 



for In ( - ) + c -> , 



So and c are constants. The matter quantities are 

4 



1 

3t2 



1 + 



In (t/r) + c 

The fluid in this case is physically self-similar. 



(5.8) 
(5.9) 

(5.10) 
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£L: The metric is plane symmetric 

ds 2 = ~dt 2 + %rS 2 dr 2 + S 2 (d6 2 + 6 2 d(j) 2 ) , (5.11) 
S r l 

where S = [ln(t/r) + c] 2 , for ln(t/r) + c — > 0, 6 and c are constants. The matter 
variables are given by 



12 4 



In ( - ) + c 



(5.12) 



Again the fluid is physically self-similar. 
• G_: In this final case the metric is plane symmetric given by: 

ds 2 = ~dt 2 + ^-dr 2 + S 2 (d6 2 + 2 d<p 2 ) , (5.13) 
S r 2 S 

where S = [ln(t/r) + c] 1//2 , for ln(t/r) + c — > 0, 6 and c are constants. 

Notice that this line element is not a solution of the system ( |4.lD -( fOD . It is just 
the asymptotic solution when ln(i/r) + c — > 0. For this reason we do not write the 
matter variables. All perfect fluid solutions approaching this singular point will tend 
to be physical self-similar. They will satisfy the energy conditions depending on the 
direction from which they are approaching this point. G_ lies exactly in the boundary 
of a region in which the energy conditions are satisfied. 

In all the solutions here presented, t, r and 9 have been rescaled in order to absorb as many 
integrating constants as possible. 



6 Special Cases 

Having completed the qualitative analysis and identified the possible asymptotic states, it 
is useful to note that in several of the invariant sets considered in the previous sections 
the system can be integrated completely so that the solutions can be written out explicitly. 
These particular sets are considered here. 
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6.1 The invariant set = 



All of the exact solutions in this case can, in fact, be determined as the system Q4.14D - Q4.15 ) 



can be integrated completely. If x% = 0, the remaining equation yields x 3 = w — x%. But 
in this case /i + p = and perfect fluid solutions are excluded. For x\ ^ we have that 

x 3 = -, (6.1) 

x 1 



and equation ( |4.15| ) becomes 

X\ = xi(wo + x\ + x\) , (6.2) 

that can be rewritten as 

{2x l + x 2 1 ) =2x l {2x l + xl + w Q ) , (6.3) 



or 



x\ — x x )' = x\(—xi + x l + wo) . (6.4) 



The case x\ = and x\ ^ implies + #i = and corresponds to the fixed points L2. 
Apart from this case the following possibilities arise. 

Case: 2x\ + x\ + w = 0. 

Notice that all the solutions correspond to the intersection of the invariants sets I2 and I4. 
The different solutions depend on the value of w . They are: 

• wq = 0: Xi = 2/x, where x = ln(t/r) + c. The metric can be written as 



ds 2 = -\dt 2 + h ^4~dr 2 + X \d6 2 + ^ 2 rf0 2 ) , (6.5) 



where b and c are constants. 

Wo = +j3 2 : Xi = — P tan(/?x/2) and 

ds 2 = -[tan(/3 X /2)] 2 dt 2 + s 2 [cos(/3x/2)] 4 ( ^dr 2 + d6 2 + sinh 2 6d(f) 2 ) , (6.6) 



so being a constant. 
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• w — —j3 2 : Two possibilities arise, Xi = /3tanh(/3x/2) or x 1 = /5coth(/5%/2). The line 
elements are 



ds 2 = -[tanh(/?x/2)]W + s^[cosh(/?x/2)] 4 I ^dr ' - dB 2 -f siir M r ) (0.7) 



and 



rfs 2 = -[coth(/?x/2)] 2 rft 2 + s 2 [sinh(/5x/2)] 4 ( ^dr 2 + d6 2 + sin 2 Odcf 2 ] , (6.8) 



respectively. 

All of these solutions satisfy the energy conditions only over some limited regions of the 
manifold. 

Case: —x x + x\ + w = 0. 

Different solutions appear again depending on wq. They correspond to the intersection of 
the invariant sets I2 and I 5 . 

• Wo = 0: xi = — 1/x, and the metric can be written as 

rfs 2 = ^ ^-dt 2 + ^rfr 2 + d9 2 + fl 2 d0 2 ) . (6.9) 

The metric is conformally flat. 

• w = +(3 2 : xi = /3tan(/3x), 

ds 2 = -[tan(/3x)] 2 rft 2 + . f Q a ... ( ^dr 2 + d9 2 + sinh 2 6dA . (6.10) 

[cos(/?x)] 2 \ r / 

• Wq — —f3 2 \ x 1 = — pta,nh(Px) or x i — ~ /3coth(/3x)- 

rfs 2 = -[tanh(/3x)] 2 rft 2 + r f? a , 12 f ^dr 2 + d# 2 + sin 2 6dA , (6.11) 

[cosh(/?x)J \ r / 



and 



ds 2 = -[coth(/^x)] 2 ^ 2 + f . *° , 19 f ^rfr 2 + d£ 2 + sin 2 6dA , (6.12) 

[smh(/5x)] 2 \ r / 



respectively. 
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In all cases b, (3 and So are constants. 

Case: 2&i + x\ + w ^ and —x\ + x\ + w Q ^ 0. 



Equations fl6.3|) and Q6.4|) can be integrated to yield 

2ii + x\ + Wq = 3aS 2 



and 



Xl + x 1 + w 



3A 1 



(6.13) 



(6.14) 



1 ' ~ u 2 S 

where a and A are arbitrary non-null constants. Comparing equations (|6.13| ) and ( |6.14j ) we 
find 



x\ = aS 2 + ^ - w • 



(6.15) 



Hence, once S is know, all the other metric functions can be calculated. ^From ( |6.1| ) we get 

(6.16) 



2$ _ 2 
6 — SC-^ . 



and also we have 



3 2 * = b 2 S 2 



(6.17) 



where b 2 = —Wok for k ^ and any arbitrary constant for k = 0. Then substituting 
x\ = S/S in equation ( |6.15|) one gets 



5 



±1 . 



(6.18) 



VaS 4 - U7 5 2 + XS 

This is the only equation that needs to be integrated. Thus, this case X2 = is completely 
solved up to quadratures. 



6.2 The case x% = (and ±3 = 0) 

This case is of particular interest since all the solutions belong to the intersection of the 
invariant sets and I§ and therefore they are physically self-similar. They represent the 
geodesic solutions for the system; i.e., these solutions have zero acceleration. Since the 
governing equations impose the condition there can be no hyperbolically sym- 

metric solutions in this case. The plane symmetric case is the special case X\ — and this 
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corresponds to vacuum solutions. Their metric is given by 

t\ 2 b 2 



ds z = -dt z + 1 + a 



: dr z + d6 z + 6 



2 i n2ji2 



(6.19) 



where a and b are constants. 



All the other solution will exhibit spherical symmetry, and the metric can be written 



as 



- 



ds 2 = -dt 2 + ^dr 2 + S 2 (d9 2 + sin 2 Odcf) . 



The governing equations reduce to 



Xi = XiX 2 , 

x 2 = x 2 (l — 3xi — x 2 ) 



and substituting x 2 from ( |6.21| ) into equation ( |6.22D , we obtain 



i\ = Sxiii + ii . 



(6.20) 

(6.21) 
(6.22) 

(6.23) 



A first integral this equation is 



(6.24) 



where n is an arbitrary constant. The special case i\ = corresponds to the curve L 2 of 
singular points. The other solutions depend on the different value of the constant n and 
they are 

• n — 0: X\ — 1/3 + 1/x; S = s (t/r) 1 / 3 x 2 / 3 , where now \ = (3/2)(ln(i/r) + c), c being 
a constant. 

• n = (3 2 : x x = 1/3 - /3tan(/3x); S = so^/rf/^cos^x)] 273 - 

• n = -/5 2 : Xl = 1/3 + /3ta,nh(Px); S = s (t/r) 1 / 3 [cosh(/3x)] 2/3 , 
or Xl = 1/3 + /3coth(/3x); S = s (t/r) 1/3 [sinh(/3x)] 2/3 . 



The solutions satisfy the energy conditions if X\ > and x% > 1/3 + 3n. 
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7 Physical self- similarity 



As was stated in the introduction, the cases of physical and geometric self-similarity are 
not necessarily equivalent when considering KSS spacetimes. Perfect fluid solutions will 
necessarily be physically self-similar if they satisfy X4 + 2x1X3 + x\ = 0, and hence lie in the 
invariant set 1^. In this case the system ( |4. 1|) -( |4.4|) reduces to the three-dimensional system 
of autonomous ODEs: 

x\ = xi(x 2 + x 3 ) , (7.1) 
x 2 = x 2 (l - 3xi - x 2 + x 3 ) , (7.2) 
£3 = x 3 (-xi +x 2 -x 3 ) . (7.3) 

Through the use of monotonic functions it can be shown that all of the asymptotic behavior 
in this class of solutions is described by solutions in one (or more) of the invariant sets 
x\ = 0, x 2 = 0, X3 = or x\ + 2x3 = 0, all of which have been previously discussed. The 
case Xi = corresponds to the vacuum case, and the exact solutions for the cases x 2 = 
and x 3 = can be found in section [TT] and respectively. 

In the case of physical self-similarity, perfect fluid solutions have a barotropic equa- 
tion of state if, in addition, they satisfy 

2x 2 + 3xi = d , (7.4) 

where c\ is an arbitrary constant. In Table f|, we summarize all the possible barotropic, 
physical self-similar solutions. 

8 Discussion 

To summarize, we have studied perfect fluid (spherically, plane and hyperbolically symmet- 
ric) space-times admitting a kinematic self-similarity of infinite type. We have restricted our 
attention to the case in which the kinematic self-similar vector field commutes with all of 
the the Killing vectors. Three different cases arise depending on the orientation of the fluid 
flow relative to the kinematic self-similar vector. The interesting general case is the 'tilted' 
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one in which the four-velocity is neither parallel nor orthogonal to the self-similar vector 
field. In this case, we have shown that the governing equations reduce to a four- dimensional 
autonomous system of ODEs. The qualitative properties of the system have been fully stud- 
ied. In particular, through an extensive use of monotonic functions, we have shown that all 
asymptotic solutions in this infinite class of kinematic self-similarity are necessarily located 
at singular points (either at finite or infinite values of the dependent variables) which are 
classified in Tables [I] and ^|. 

Most of these singular points are saddle points in the full phase space, although there 
are invariant regions in which they do act as sinks or sources, thereby acting as attractors 
(or repellors) for classes of solutions. The only global sinks and sources are located on the 
infinite boundary, summarized in Table |5|. Hence, in general solutions asymptote to one of 
those represented by the points A + , B + or D_ in the past and one of those represented by 
the points A_, B_ or D + in the future. 

The physical solutions described by these singular points are given in the cases in 
which the weak and dominant energy conditions are satisfied. The class of solutions which 
are also physically self-similar are again important in this analysis. We show that in all cases 
in which the energy conditions are satisfied the asymptotic behavior is necessarily physical 
self-similar and the space-time is plane or spherically symmetric. This result coincides with 
the results of Benoit and Coley ITOj , which studied the case of spherical symmetry with 



finite kinematic self-similarity. This again shows the relevance of the physical self-similar 
models. 

In some special cases, e.g., the invariant set £2 = 0, and the geodesic case, corre- 
sponding to the case £3 = 0, the four- dimensional autonomous system of ODEs can be 
integrated completely. All the exact solutions have been found in these cases (see Sections 
6 and 7). In the geodesic case, the solutions are again physically self-similar. These exact 
solutions serve as illustrations of the more general qualitative results previously discussed. 
Finally, we have also found all of the physical self-similar solutions that admit a barotropic 
equation of state. The results are summarized in Table 4. 
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Table 1: Finite Singular Points for equations ( (HQ - flOp . The local analysis for each 
singular point will be discussed in the subsequent sections, according to their classification. 
Qi,2,3 are isolated singular points and Li and L 2 are curves of non-isolated singular points 
(one-dimension equilibrium sets). Note that the curves L 1 and L 2 intersect at the point 
(0, 0, 0, 0), which is a saddle-node. (3 and a are constants. 



(^1,^2,^3,^4) 


Eigenvalue - Eigenvector Pairs 


Classification 


Qi (0,1,0,0) 


1 (0,1,2,0) 

1 (1,0,3,0) 
-1 (0,1,0,0) 

2 (0,1,3,3) 


Saddle 


Q 2 (1,-1,1,0) 


-3 (5,-6,-9,0) 
-2 (-4,5,3,7) 

1 + V2 (1,^,1,0) 
1-V2 (1,-^,1,0) 


Saddle 


Q3 ill \ 1 0) 


3/4 (1,0,3,0) 
1/4 (-2,1,-3,2) 
-(l->/7i)/8 (-2,-^,1,0) 
-(l + V7i)/8 (-2,^,1,0) 


Saddle 


Li (0,0,/3,/? 2 ) 


(3 (3,0,1,0) 
1 + /3 (0,l + /?,/3,2/? 2 ) 
-2/5 (0,0,1,0) 

(0,0,1,2/3) 


Curve of Saddle Points 


L 2 (a, 0,0, -a 2 ) 


(-1,0,1,0) 
2a (1,0,2,0) 
l-3a (a, -2a, 0, -3a) 
(1,0,0, -2a) 


Curve of Saddle Points 
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Table 2: Classification of the Infinite Singular Points of the Poincare transformed 
system ( [4.1| )-( f4~4|) 1. The exact solutions for each of these points are described in section 5.2. 



(Xi, X2, X3, X4) 


Eigenvalue - Eigenvector Pairs 


Classification 


A± (0, ±1,0,0) 


±2 (1,0,0,0) 
±2 (0,1,0,0) 
±2 (0,0,1,0) 
±3 (0,0,0,1) 


A + : Source 
A : Sink 


B ± (0,0, ±1,0) 


±2 (1,0,0,0) 
±2 (0,1,0,0) 
±2 (0,0,1,0) 

±1 (0,0,0,1) 


B + : Source 
B_ : Sink 


C± (0,0,0,±1) 


(1,0,0,0) 
(0,1,0,0) 
(0,0,1,0) 
(0,0,0,1) 


Saddle- Node 


D ± ^(±1,0, ±1,0) 


=F<\/2 (1,0,1,0) 
T^y (0,0,0,1) 
(-1,0,1,0) 

(0,1,0,0) 


D+ : Sink 
D_ : Source 


E± ^(T2,0,±l,0) 


±^fi (0,1,0,0) 

T^f (1,0,2,0) 
(0,0,0,1) 
T 2 ~f (-2,0,1,0) 


Saddle 


F ± ^(0,±1,±1,0) 


TV2 (0,-1,1,0) 
±y/2 (2,-1,1,0) 
±v / 2 (0,0,0,1) 

(0,1,1,0) 


Saddle 


G ± #(T2,±3,±1,0) 


±43 (0,0,0,1) 
T ±M (2,-3,-1,0) 

T ^m (1,0,2,0) 

T 3yli (0,1,-3,0) 


Saddle 


H ± ^(±1,^3/2, ±1,0) 


(134,-201,257,0) 
±# -3, 2,0) 
±^P (63,118,63,0) 

T zhm (0,0,0,1) 


Saddle 



Table 3: Invariant regions of the space (0:1,0:2,0:3) for the system ( 4.19|) -( 4.21|) with 



corresponding monotonic functions. Note that the sets Ui 3 — U\§ are not included here as 
they are two-dimensional invariant sets and their complete dynamics have been summarized 



s Label 


Definition of Region 


Monotonic 
Function 




{(x 1 ,x 2 ,x 3 )\x 1 ,x 2 ,x 3 > 0} 


Xi strictly increasing 


u 2 


{(xx,x 2 ,x 3 )\xi,x 2 > 0,xi + 2x 3 > 0} 


Xi + 2x 3 strictly increasing 


u 3 


{(xx,x 2 ,x 3 ) \xx,x 2 > 0,xi + 2x 3 < 0} 


x\ + 2x 3 strictly decreasing 


u 4 


{(0:1,^2,^3)1^1 < 0,a;2 > 0,xi - x 3 < 0} 


x\ — x 3 strictly decreasing 


U 3 


{(xi, x 2 ,x 3 ) \xi <0,x 2 > 0,xi - x 3 > 0} 


x\ — x 3 strictly increasing 


u 6 


{(x t , 0:2,0:3)12:1 < 0,^2,0:3 > 0} 


x\ strictly decreasing 


U 7 


{(0:1, 0:2,0:3) \xi > 0,x 2 < 0,xi - x 3 < 0} 


X\ — x 3 strictly increasing 


Us 


{(xi,x 2l x 3 )\x 1 > 0,x 2 < 0,xi - x 3 > 0} 


X\ — x 3 strictly decreasing 


u 9 


{(x 1 ,x 2 ,x 3 )|x 1 > 0,x 2 < 0,x 3 < 0} 


x\ — x 3 strictly decreasing 


u l0 


{(x 1 ,x 2 ,x 3 )|a; 1 ,X2,X3 < 0} 


X\ strictly increasing 


U u 


{(xi, x 2 , x 3 )\xx, x 2 < 0,xi + 2x 3 < 0} 


X\ + 2x 3 strictly increasing 


u 12 


{(xi, x 2 , x 3 ) Xx,x 2 < 0, xi + 2x 3 > 0} 


X\ + 2x 3 strictly decreasing 
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Table 4: Infinite kinematic self-similar perfect fluid solutions admitting a 
barotropic equation of state. The weak and dominant energy conditions are always 
satisfied, except in the last case in which c\ > 1. 



Case 


Solution 


Remarks 


ci = 


X\ = 0, X2 = 0, X3 = [ln(t/r) + c] -1 , x 4 = 


vacuum 


Cl = 1 


xi = 1/4, x 2 = 1/8, x 3 = -1/8, x 4 = 


Q%: stiff-matter 


Cl = 1 


xi = 1/3, x 2 = 0, X3 = 0, X4 = —1/9 


L 2 : stiff- matter 


ci = 2 


xi = 0, X2 = 1, X3 = 0, X4 = 


Q\\ vacuum 


ci = 2 


xi = 2/3, x 2 = 0, x 3 = 0, x 4 = -4/9 


L 2 : dust 


ci = 2 


X3 = 0, x 4 = — Xi, X2 = 1 — 3xi/2 
see case x 3 = 0, with /3 = 1/3, n — — 1/9 


dust 


ci ^0,1,2 


xi = ci/3, x 2 = 0, x 3 = 0, x 4 = -cj/9 
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Table 5: Summary of the global asymptotic behaviour. 



Singular Point 


Nature of the Singular Point 


Physical Characteristics of Solutions 


A + 


Source: Global attractor to the past 


Geodesic, vacuum solution; 
physically self-similar 


A 


Sink: Global attractor to the future 


Geodesic, vacuum solution; 
physically self-similar 


B+ 


Source: Global attractor to the past 


Vacuum solution; 
physically self-similar 


B_ 


Sink: Global attractor to the future 


Vacuum solution; 
physically self-similar 


D+ 


Sink: Global attractor to the future 


Energy conditions not satisfied 


D 


Source: Global attractor to the past 


Energy conditions not satisfied 
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Figure 1: Phase portrait of the Poincare transformation of system (|4.14 )- ( 14.151) : The case 
Wq < corresponds to spherically symmetric solutions, w$ = to plane symmetric solutions, 
and Wq > to hyperbolically symmetric solutions. The vertical direction represents the 
parameter space wq G 3?, where wq = is the bifurcation value for the system. 
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Figure 2: Dynamics on the infinite boundary for the plane symmetric solutions: for 

Poincare transformed space. Top hemisphere; i.e. X 2 > 
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Figure 3: Dynamics on the infinite boundary for the plane symmetric solutions: for 

Poincare transformed system. Bottom hemisphere; i.e. X 2 < 
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Figure 5: Phase portrait of the Poincare transformed system ( f4.24|) - (|4.25|) : plane sym- 
metric solutions restricted to the invariant set x± + 2^3 = 0. The phase space is X 2 vs. 

x 3 . 
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